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^ ■ Abstract 

It was argued many years ago that translational symmetry breaking due to the appearance 
T^ \ of spin-Peierls ordering (or bond-charge stripe order) is a fundamental property of the quantum 

paramagnetic states of a large class of square lattice antiferromagnets. Recently, such states were 
^ ■ shown to be a convenient point of departure for studying translational symmetry breaking in 

doped antiferromagnets: these results are briefly reviewed here with an emphasis on experimental 

implications. In the presence of stronger frustration, it was also argued that the insulating 

antiferromagnet can undergo a transition to a deconfined state with no lattice symmetry breaking. 
^^ ' This transition is described by a fully-frustrated Ising model in a transverse field: details of this 

Q ■ earlier derivation of the Ising model are provided here — this is motivated by the reappearance 

of the same Ising model in a recent study of the competition between antiferromagnetism and 

d-wave superconductivity by Senthil and Fisher. 
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I. INTRODUCTION 

Following the fundamental neutron scattering observations of Tranquada 1], there has 
been a vigorous resurgence J2, |3[ in the study of states of two-dimensional strongly correlated 
systems which break various types of translational symmetries (to be precise: we will consider 
translational symmetry to be broken when an observable invariant under spin rotations (and 
which does not change the net charge of the system) does not remain invariant under the 
space-group of the lattice). In an early discussion of translational symmetry breaking in two- 
dimensional antiferromagnetsj^, |5| it was argued that any quantum paramagnet accessed 
by a continuous quantum transition from a Neel state must have (for S = 1/2 spins) (i) 
broken translational symmetry due to spin-Peierls (or bond-charge density) ordering, and 
(m) confinement of spin-1/2 excitations resulting in integer spin quasiparticle modes. Other 
early work[6|,|7| used Hartree-Fock theory to study the formation of site-charge density order 
in lightly doped antiferromagnets. In this paper, we shall review recent work[8|, |9| in which 
these two apparently divergent physical effects were studied in a combined footing — we 
shall argue that such an analysis leads to considerable physical insight and has significant 
experimental implications. 

This paper is organized as follows. In Section |n] we will consider translational symmetry 
breaking in the quantum paramagnetic phases of insulating antiferromagnets — this will be 
done using lattice 'height' and Ising models which are obtained after a duality transforma- 
tion. In Section UTTl we will turn to the recent results^, iQi] on doped antiferromagnets, and 
their implication for experiments. 



II. INSULATING ANTIFERROMAGNETS 

A convenient and physically transparent framework for discussing the physics of the low- 
lying singlet excitations of a quantum paramagnet is the "quantum dimer model" [iCf]. At 
first glance, the formulation of the dimer model does appear rather arbitrary, with a number 
of ad hoc assumptions which obscure the connection to the Heisenberg antiferromagnet — 
instead, the dimer model is best viewed as a caricature which captures some essential pieces 
of the physics. We shall mainly use the dimer model as an intuitive intermediate step 
towards introducing dual 'height' and Ising models. These dual models have also been 
derived by a number of other methods, including systematic semiclassical and large N 
expansions!^, |5|, and these provide the more formal justifications for the arguments here. 
We will describe duality transformations on an effective model of resonating singlet bonds 
which show that unfrustrated antiferromagnets are described by an effective height model 
in 2+1 dimensions[5. 12lll3l|: in the presence of frustration, the height model gets modified 
to a fully frustrated, two-dimensional Ising model in a transverse field 1J| . All of the results 



in this section were derived and published some time ago|5l.ll2Lll3lll4l| ; the frustrated Ising 
model was derived by one of us'lS' and the main results were outlined in Ref.[14j , but de- 
tails have not been previously published. The presentation of the details here was motivated 
by the reappearance of the same Ising model in a recent study of the competition between 
antiferromagnetism and d-wave superconductivity by Senthil and Fisher 16|. We also note 



another recent study of lattice gauge theory models of two-dimensional antiferromagnets by 
Nagaosa and Lee[l7|; they examine theories related to those considered in Ref(l3| and in the 
present paper, but their actions do not contain the crucial Berry phase terms (as in (j2.16|) . 
(j2.17|) . and (|2.25|) ) associated with the non-zero net matter density on each site even in the 
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FIG. 1: Dimer coverings of the square lattice, each representing a distinct state in the singlet Hilbert 
space. The two states in (a) can resonate with each other as can the two in (b). The columnar 
state in (a) has the niaxinium number of resonating partners, and this drives the translational 
symmetry breaking. 



undoped antiferromagnet — it is these Berry phase terms which are responsible for frustra- 
tion in the dual height or Ising spin representation, and we believe their incorporation is 
essential for a proper description of the physics. 

The central assertion of the dimer model is that it is appropriate identify the Hilbert 
space of the low-lying singlet excitations of the quantum paramagnet on a square lattice 
with a set of dimer coverings lOj. Each dimer covering corresponds to different singlet 
state, and off-diagonal matrix elements between such states therefore lead to a "resonance" 
between the dimers. Such resonance terms must preserve the constraint of one and exactly 
one dimer terminating at every site, and this severely restricts their form. Two possible 
resonating terms are illustrated in Fig^ One might now anticipate that the ground state 
of such a resonating-valence-bond Hamiltonian is a liquid-like superposition of all dimer 
coverings, but this turns out not to be the case. The ground state in fact breaks lattice 
translational symmetryjj, |5| e.g. for some range of parameters it looks like the columnar 
ordered state at the top of Fig ^ There is a simple intuitive reason for this: the columnar 
state has the maximum number of plaquettes about which the dimers are able to resonate 
and still preserve the local constraints. This resonance lowers its energy, and this effect is 
strong enough that the system always picks out one of the four columnar states and retains 
memory of its orientation. The ordering is therefore due to a quantum "or der-from- disorder" 
effect. 

Let us now make these arguments more precise. We will divide our discussion into two 
subsections. Section 111 Al will consider the model with dimers connecting nearest neighbors 
on the square lattice. More generally, the results here will apply to bipartite lattices in which 
the dimers only connect sites on opposite sublattices. We will show that such models here 
are described by effective height models in 2+1 dimensionsp, Uj. Section fllBI will consider 



models in which the dimers are allowed to connect sites on the same sublattice — the effective 
model for these will be the fully frustrated Ising model in a transverse field IJ . 



A. Bipartite dimer models 

We will only discuss the case of the nearest-neighbor, square lattice dimer model here. 
We identify the dimers by site, i, of the square lattice on their lower or left end. Let rjiEia 
be the number operator for the dimer on site i oriented in the a direction (a = x,y); rji 
indicates the sublattice of site i, and equals +1 on one sublattice and —1 on the other. 
We have so far considered only a S = 1/2 antiferromagnet, but for general S there will be 
exactly 25" dimers emerging from every site — this constraint translates to 

A„E,, = 2%, (2.1) 

where A^^ is the discrete lattice derivative in the a direction. The factors of rji where 
introduced so that the constraint would have the Gauss-law form in (j2.1|) . and the operator 
Eia is seen to be the analog of the electric fieldjl^]- We also introduce an angular phase 
variable, Ajq, (the analog of a compact f/(l) gauge field), on every link which is canonically 
conjugate to Eia. 

[Aia, Ejp] = i5ij5a/3, (2.2) 

where it should be clear from the context when we mean i = \/— T, and when i is a site 
label. 

We can now write down the Hamiltonian of the dimer model: 

H, = ^J2eL-K,Y: cos(e,^A,i,^). (2.3) 

i,a i 

The first term, proportional to Ki is only non-trivial when 25* > 1, and it ensures that the 
density of dimers is as uniform as possible. It follows from the commutation relations ()2.2|) 
that the second term, proportional to K2, flips dimers around a plaquette and so induces 
the 'resonance' shown in Fig^ 

We will now write down a path integral representation of the partition function of Hd by 
following a standard route. We insert complete sets of Ei^ eigenstates at small imaginary 
time intervals Ar. The matrix elements of the cosine term in Hd are evaluated by replacing 
it with the Villain periodic Gaussian form (this is the only 'approximation' made in our 
duality mappings - all other transformations below are exact): 

exp {K2/^t cos{ea(i /^aAii3)) ^ V exp - " h iBatap/^aAip ■ (2.4) 

V ^ B, \ 2K2AT J 

Here Ba is an integer-valued fleld on the sites, a, of the dual lattice; there is an obvious 
relationship between the location of the direct lattice site, i, and the dual lattice site a, but 
we will not specify it explicitly to avoid cluttering up the notation (we will consistently use 
the labels i,j... for sites on the direct lattice, and the labels a,b. . . for sites on the dual 
lattice). In the present case, a resides at the center of plaquette around which the 'flux' 
eQ/jAoAj^ resonates the dimers. It is also convenient to introduce a three-vector notation in 
space time: we deflne the integer-valued 'electromagnetic flux' vector F^^ = {Eiy, —Ei^, —Ba) 



• ,••,•• ••••• 

r^ X X^^ X 1/4 1/4 

• .••.•• ••••• 

>^-' X >^"'' X 3/4 1/2 3/4 1/2 

•k'^ X >^+^ X 1/4 1/4 

• •••• ••••• 

(a) Xa, (b) ya 

+ 1/8 -1/8 +1/8 -1/8 

XXX 
-1/8 +1/8 -1/8 +1/8 

XXX 

+1/8 -1/8 +1/8 -1/8 
(c) Zi^ 



FIG. 2: The values of the only non-zero components of the fixed field Xa^, ya, and Zi^. The 
circles (crosses) are the sites of the direct (dual) lattice. In (c), only the fj, = t component of Zj^ 
is non-zero and its values are shown. 

on the dual lattice sites, where the index /x = (x, y, r) (we will consistently use the labels 
a, /3 ... to represent spatial components only, while /x, z/, A . . . will represent three-dimensional 
spacetime components) . Here Eia refer to the integer eigenvalues of the operator Eia which 
are summed over in each time step. Now the partition function of H^ can be written in the 
compact form 

Zf= Y: exp (-^T. fA n S {e,,xKFax - 2%5,,) (2.5) 

Here the sum is over the integer-valued field F^^ which resides on the sites of the dual 
cubic lattice in spacetime; the delta function constraint imposes 'Gauss's law' ()2.1|1 and the 
'Maxwell equation' between E and B which follows from the matrix elements of ()2.4|1 . The 
time spacing At has been chosen so that the coupling e^ = KiAt = 1/{K2At). 

We now solve the constraint in ()2.5p by writing F^^ as the sum of a particular solution 
and the general solution of the homogeneous equation: 

F„^ = A^Na + 2SXa^,. (2.6) 

Here A^^ is a fluctuating integer-valued field on the dual lattice sites, while Xa^ is a fixed 
field independent of r satisfying 

ef^^x^u^ax = r]i6f,r. (2.7) 

A convenient choice is to take Xax = 0, Xar = 0, and Xay as shown in Fig|2K, taking the 
values ±1 on every second column of sites and zero otherwise. For future manipulations, it 
is convenient to split Xafj, into curl-free and divergence-free parts by writing 

^afM = ^fj.ya + ^iiuX^u^iX, (2.8) 



where again 3^^ and Zi^ are fixed fields independent of r and their values are shown in 
Fig|2b,c; ya takes the values 0, 1/4, 1/2, 3/4 on the four dual sublattices, while Zj^ = 6^rVi/^- 
We have now assembled all the ingredients necessary for the final representations of 
the partition function. There are three different formulations, and each present a valuable 
physical perspective and can be used for further quantitative analysis. We will present them 
in turn in the following subsections. 

1. Height model 

We insert fl2.6|) into ()2.5|) and obtain the height model partition function 

Zh=Y. exp ( -^ E i\^Na + 2SXa^)A . (2.9) 

For the purposes of this section, it is convenient to insert the decomposition ()2.8|1 in (J2 



ror tne purposes oi 
and obtain iii. El 

Z'u= Eexpf-^E(A.^a)')' (2-10) 

where the 'heights' Ha are defined by 

if, = iV, + 25X. (2.11) 

Notice that Z^ has dropped out-this is a general property of all representations of bipartite 
dimer models, but the non-bipartite models in Section FlI Bl will depend upon Za- We can 
view Ha as the heights of a 3-dimensional interface which are restricted to take values equal 
(for S = 1/2) to those in Fig|2|D plus arbitrary integers. It is instructive to note here, for 
S = 1/2 , a direct connection between these heights and the original dimer model we began 
with[l8|. For large e^, the partition function Z'jj will be dominated by configurations in 
which neighboring heights are as close to each other as possible; the manifold of minimum 
action is defined by heights which satisfy 

\Ha — Hh\ < 1 for every pair of nearest neighbors a, b, (2.12) 

If we consider two height configurations to be equivalent if they differ only by a constant 
integer shift Ha — > Ha+p, then the set of all such equivalence classes is precisely equivalent to 
the set of dimer coverings. The identification is illustrated in FigEl neighboring heights differ 
either by large (3/4) or small (1/4) steps, and we place a dimer on every link intersecting a 
large step. 

2. sine- Gordon model 

Promote the integer valued field Na to a real-valued field (fa by the Poisson summation 
formula. After shifting the real field by 'fa^'^a — ya we obtain from (J2.5I2.6|2.8|) 

/oo / p2 \ 

d^a E exp --E (A^^a)' + 2mY,Pa {^a " 3^a) • (2.13) 

'°° {Pa} \ "^ a,fi a J 



{Pa} 

If we now truncate the summation over Pa by adding a fugacity e^^^^^^" , then 

-iPa 



^ ^2niPa{^a-ya)+iiny)PS ^ ^2ycos(2niva-ya))^ and ZsG is seen to be a version of the sine-Gordon 



field theory in 2+1 dimensions!^, |5[. 
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FIG. 3: A set of values of Ha satisfying H2.12() and the corresponding dimer covering. 

3. Compact QED 

Define M^^ = A^N^, so M^^ is an integer- valued field satisfying the constraints 

e^,xA,Max = 0. (2.14) 

We impose these constraints by a real Lagrange multiplier field Ai^ by inserting factors of 

— ^ exp {iAi^e^,xA,Max) (2.15) 

-IT Zn 

in Zh in ()2.9|1 . The sum over the Mq^ can now be performed using the Poisson summation 
formula and we obtain after using ()2.7|1 



ZeQED= E n/''^^exp|-^^(e^,AA,A,A-27rg„^)V^255:r/,A,J, (2.16) 



where the Q^^ extend over all integers. Z^qed is easily seen to be the Villain form of the 
partition function of a compact U{1) gauge theory in 2+1 dimensions. There is an additional 
Berry phase carried by the world lines of the background charge density of 25* particles on 
each site. An alternative, and exactly equivalent, form of Zcqed can be obtained by applying 
the above transformation to Z'^j (in ()2.10|) ) rather than to Zh- Then, exactly the same steps 
lead to 



/•TT f]A. / 1 

Z'cQED= E ny_^^exp -— ^(e^,,A,A,,-27rg,^)Vz47r5^3;,A^Q, 



r-^ dAi„ ( 1 

j^[e^,^xAyAix-2T\Qai,y ^^^'nS^yaA^S^a^l 

(2.17) 
The instanton number of the gauge theory is A^Qa^, and to the Berry phase is now carried 
by the instantons: each instanton has a Berry phase of e^^'"^^°-. It is a remarkable fact 
that precisely the same Berry phase was obtained in very different semiclassical and large A^ 
analysis of bipartite antiferromagnets^, JJ:]. It is also worth reiterating that Zcqed (in which 
the Berry phase is carried by the world hues of the 25* particles on each site) and Z^p ^n (in 
which the Berry phase is associated with the instantons) are exactly equivalent J13l| — they 
merely do the book-keeping of the total Berry phase somewhat differently. 

The properties of the theories Zo-, Z_^n, and Z'^q^j^ are reasonably well understood and 
have been described elsewhere|5l. Il3l. Ilq . The height model is in a smooth phase phase for 
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FIG. 4: Another possible ground state of the dinier niodeh the plaquette state. 

all values of e^ (there is no roughening transition), while the compact U{1) gauge theory 
is always confining. In terms of the underlying dimer model, this means that translational 
symmetry is spontaneously broken and the dimers crystallize in a bond-charge ordered state. 
A likely configuration is the columnar state in Fig ^, but it is not ruled out that one of 
these theories could exhibit a transition to some other confining state in which the details of 
the lattice symmetry breaking are different: another state^] which occurs in the parameter 
space of ZsG is the plaquette state of Fig HI 

B. Non-bipartite dimer models 

We will now relax the constraint that the dimers must connect sites on opposite sublat- 
tices. We will show that such a step reduces the height model to a 2+1 dimensional fully 
frustrated Ising model in a transverse field, following the derivation (iq outlined in Refill] . 

An example of a dimer state connecting sites on the same sublattice is shown in Fig |^. 
Notice that a pair of such 'diagonal' dimers can resonate with a pair of dimers connecting 
opposite sublattices. However models which include terms like those in FigEt^ do not appear 
to be amenable to the exact duality transformations we have considered in Section IIIAl So 
we have two alternatives for further analysis: we could be satisfied by approximate duality 
transformations on models which include the resonance in FigEfe^, or we could modify the 
underlying dimer model to allow duality to proceed smoothly. Both alternatives lead to 
essentially the same final result, but we will choose the latter for presentation here as it 
more clearly isolates the approximation made. The modification of the dimer model is 
aided by the physical picture developed in the large- A^ theories of the effects of frustration 
presented in Refs J2ll . |22| . The key property of the diagonal dimers in Fig \^ is that they 
connect sites on the same sublattice, and such dimers carry net charge ±2 under the compact 
QED representation of Section III A 31 This will be shown explicitly in more detail below, 
but can be seen quickly as follows. From the representation ()2.16|) we see that the model 
behaves as if there are background charges of +2S on one sublattice, and charges of —25" 
on the other sublattice. A dimer connecting sites on opposite sublattices absorbs (in the 
Gauss Law constraint (J2.H) ) charges of ±1 on its two ends; so such a dimer has net charge 
0. However we do need a vector potential for a gauge invariant connection between the 
charges on the two ends, and this is why the number operator of the dimer ~ e*"^*" (see 
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FIG. 5: (a) Dimers which connect sites on the same sublattice resonating with dimers connect- 
ing sites on opposite sublattices. (b) The same sublattice dimers have now been collapsed onto 
monomers representing singlet bonds between electrons on the same site. We argue that the defor- 
mation from (a) to (b) should have no significant effect on the long distance properties: it is crucial 
to keep track of which sublattice the two ends of a singlet bond reside, and local deformations 
which preserve this information will not modify the long distance physics. The diagonal dimers in 
(a) and the monomers in (b) both carry net charge ±2, and this is the key factor which controls 
their effect on the long distance physics. 

()2.2|) ). Let us now apply the same argument to the diagonal dimers of Fi^S^. These connect 
sites on the same sublattice and therefore carry net charge ±2; however, we also need a 
vector potential to connect such spatially separated charges, and so the number operator of 
a diagonal operator must also contain the necessary path integral of Ai^. It is this combined 
role of diagonal dimers that precludes a straightforward duality transformation, but it also 
suggests a simplifying modification: the important thing is that the dimer carry charge ±2 
while preserving the fact that each site has net charge ±25"; so we collapse the two ends of the 
diagonal dimers onto the same site, as in FigEjD. The resulting monomer trivially satisfies 
the requirement of representing a singlet bond between electrons on the same sublattice 
and thus carries charge ±2; moreover, no vector potential is required as the charges on its 
two ends are on the same spatial location. Further, such monomers can resonate, like the 
diagonal dimers, between dimers connecting sites on opposite sublattices, as shown in Fig El 
Of course such a model requires that the number of singlet bonds allowed to emerge from 
a site be at least two i.e. 2S > 2. Nevertheless, we maintain that generic properties of 
the model studied here apply also to S* = 1/2 models with configurations such as those in 
Fig Eli. 

Let us now write down the Hamiltonian of the monomer-dimer model containing the 
resonance in FigO). Let rjifii be the monomer number operator on site i. Then the constraint 
that there are exactly 25" singlet bonds emerging from site i modifies the 'Gauss Law' 
constraint ()2.1|) to 

A^Ei^ + 2hi = 2Sr]i. (2.18) 

This equation shows that in addition to the background 'charge' of 25", there isJ2l|, |22|, as 
claimed above, a fluctuating charge density of a 'Higgs field' of monomers of charge 2. We 



also introduce a angular phase variable 0j canonically conjugate to fii. 

[(puhj] = i6ij (2.19) 

We can easily now extend the dimer Hamiltonian H^ in (j2.3|) to include the resonance in 
FigElD, induced by the following cosine term proportional to K^: 

Hd = ^T.^ia + ^J2^^i - K^T. cos(e«;3A„ii^) -K,J2 cos(A,0, - 2Ai^). (2.20) 

i,a i i i.a 

For convenience we have also added an additional term, proportional to K^, which controls 
the density of monomer fluctuations. ^ 

The remainder of this subsection describes the analysis of the Hamiltonian Hd, with the 
quantum dynamics defined by the commutation relations ()2.2j) and ()2.19|) . while preserving 
the local constraint (j2.18|) . The procedure will closely parallel that followed in Section 111 Al 
We insert complete sets of Eia and Ui eigenstates at intervals of imaginary time At in 
the partition function. As in ()2.4j) . we replace the exponential of the new cosine term 
proportional to 7^4 by a Villain form but with a decoupling field ki^. After evaluating the 
matrix elements between the complete sets of states, defining the three- vector integer- valued 
'current' jj^ = {ki^, kiy, —fii) we find that Zp in (j2.5j) is replaced by 

Zf= E expl-^Y.K-^T.J^^]lH^,uxKFaX + 2J.,-2Sv^S,r)l[^{A,J,,), 

(2.21) 

where, for simplicity, we have chosen the coupling g = K^Ar = l/{K4Ar). The sum in 
()2.21|) is over integer valued fields Fa^, on the dual lattice, and ji^ on the direct lattice. The 
solution of the constraints analogous to ()2.6|) is 

Fa^ = A^N, - 2a,^ + 2SXa^, (2.22) 

where Oq^ is an integer-valued 'gauge' field on the links of the dual lattice, while, as before, 
Na is an integer-valued height field on the sites of the dual lattice. Inserting (|2.22p into 
()2.21|) we obtain the generalization of the height model ()2.9|) 

Zh= E exp ( -^ E i^f^Na - 2a,^ + 2SXa^,f - | ^ {e^uxA^aaxf] • (2.23) 

We will refer to this partition function as the 'gauged height model'; the form ()2.23j) is the 
most convenient for subsequent analysis and is one of the key expressions of this paper. We 
also note an alternative form of Zh, analogous to ()2.1()|1 . obtained by inserting ()2.8|1 into 

Z'h= E exp I -^ E (^M^a - 2aaM)' " f E UuX^'^<^aX - '^Z,,] I , (2.24) 

where the heights Ha are defined in ()2.1H1 . As promised earlier, Zi^ does not drop out of 
the models of this subsection. 

The main remaining task is to describe the physical properties of the gauged height 
model ()2.23|) : we will embark on this in Section FlI B 21 However, before we do so, we will, 
for completeness, present the analog of the compact QED models of Section III A 31 

10 



1. Compact QED 

We apply the method of Section III A 31 to Zh- We define integer- valued fields M^^ = 
A^A^a - 2aa^ and 6^^ = e^^x^u^ax which satisfy the constraints e^^x/^^Max + 26^^ = and 
^fibifj, = 0. Imposing these constraints by continuous real fields Aj^ and 0j respectively, and 
performing the summation over M^^ and 6j^ by the Poisson summation formula we obtain 



^cQED= E n/_ ^n/ ^e^P -^E(eM^AA^Aa-27^Q„^)^- 

^ J2 (A^^iM - 2^.M - ^^Pif^? + ^2^ E ^i^- I ' (2-25) 



1,11 



where Qafj. and pi^ extend over all the integers. This is the Villain form of the action of 
compact QED coupled to a charge 2 Higgs scalar |^,|2^. The Berry phase term, attributable 
to the background charge density of 2Srii on each site, is the same as that in Zcqed in ()2.16|) . 
If, as in Section I11A3| we had applied the transformation of this subsection to Z'^ rather 
than Zh, then the Berry phases would have been attached to the monopoles and vortices of 
the compact QED theory 1J| ; this is straightforward to do and we will not explicitly present 



the results here — the final form is, of course, exactly equivalent to ()2.25p . ()2.24|) or ()2.23p . 
Also note that the g = oo limit of Zcqed in ()2.25p is exactly equal to Zcqed in fl2.16|) . 

2. Gauged height models 

We now return to the main objective of this section: determination of the phase diagram 
of the gauged height model Zh in fj2.23p . We will restrict our attention to S" = 1/2; results 
for other half-odd-integer values of S are similar, and the generalization to integer S is 
straightforward. 

As in the work by Fradkin and Shenker [2j|, it is useful to look at various limiting values 
of the couplings in Zh' from these results we have constructed the phase diagram shown in 
Fig El 
(z) g = oo: 

The field a^^ must be 'pure gauge' along this line, with a^^ = A^Cq for some integer-valued 
Ca- This can be absorbed in Na by Na —<■ Na + 2ca, and so Zh reduces to the pure height 
model Zh considered in (|2.9|) in Section 111 Al So the heights are in a smooth phase, the 
gauge theory is confining, and there is translational symmetry breaking due to the presence 
of bond-charge-density order. 
(ii) e^ = 0: 

Now the height field Na drops out, and the theory is simply that of the integer- valued gauge 
field Qafj. with a Maxwell action. This theory is exactly equivalent to the pure XY model 
in 2+1 dimensions, as can be seen by setting Ai^ = Qa^ = in ()2.25|) . So there is a phase 
transition in the D=3 XY universality class at a critical g = gc, the g < g^ {g > gc) phase 
maps onto the low (high) temperature phase of the XY model. 
(m) g = 0: 

This is the most interesting limit, where new physics emerges. With no Maxwell term for 
the gauge field Oa/^, the sum over a^^ can be performed independently on each link. From 

11 



this it is evident that only the parity of the Na is relevant-the partition function depends 
only upon whether a given Na is even or odd. So we can introduce an Ising spin variable, 
aa on each site of the dual lattice by defining 

aa = 2 [Naimod 2)] - 1. (2.26) 

In terms of the aa, Z^ reduces to the Ising partition function in three dimensions [ij] 



Zl = Yl ^^P ~ H JabCTaO-b 

{<Ta=±l} V <a,b> / 



(2.27) 



The sum a, h is over nearest neighbor pairs on the dual cubic lattice. The exchange constants 
J ah all satisfy \Jab\ = -^(e^) where J(e^) is a monotonically increasing function defined by 



.2J{e2) ^ 



.n=— oo 



-2n2e2 



/ 



.n=— oo 



-{2n+l)2e2/2 



(2.28) 



The presence of the Xa^ in Zj:^ introduces a key variation in the signs of the Jab- it is not 
difficult to see that these signs must be chosen so that each plaquette in the x-y plane 
is frustrated 1.4|. In the language of 2+1 dimensional quantum models, this is the fully 
frustrated Ising model in a transverse field J2J]. Such a model has been studied in earlier 
works J^,^,^], and we summarize the main results. There is a phase transition at 14. 26| 
e^ = e^ where J(e^) ~ 0.35, such that for e^ > e^ both the Ising symmetry iX'^,^ ^ 0) and 
the translational symmetry are simultaneously broken. The broken translational symmetry 
is due to development of columnar spin Peierls order of the type shown in Fig ^. In the 
language of the gauge theory, Zcqed, this phase is confining. For e^ < e^, the gauge theory 
is in the deconfined 'Higgs' phase, and Ising and translational symmetries are restored. This 
mechanism of deconfinement of spinous in antiferromagnets by condensation of a charge 2 
Higgs field was first proposed in Ref 2l| . The condensation of the Higgs field means that 
there are strong fiuctuations in the conjugate number operator i.e. in the number of dimers 
carrying charges ±2 — these can fiuctuate into dimers of net charge 0, without violating the 
number constraint which fixes the total number of singlet bonds from each site at ±25", as 
shown in Fig \^ According to the arguments of Refs [ij, |25| the phase transition in the fully 
frustrated Ising model at e^ = e^ is also in the XY universality class, but now the e^ > e^ 
(e^ < e^) region maps onto the low (high) temperature phase of the XY model. Finally 
we note that the simulations of Grest J26| appear to observe a second phase transition at 
e^ = e|^ where J{e\^ ~ 1.8: for e^ > e\^ the pattern of the translational symmetry breaking 
appears to change into the plaquette phase of Fig |3] — however the gauge theory remains in 
the confining phase. 



,iv) e 



cxd: 



The physics here is very similar to the large e^ region considered in [Hi) at g = — 
translational and Ising symmetries are broken and the gauge theory is confining. 

All of the above results have been combined in Fig El which is one of the central results 
of this paper. There is a single phase boundary connecting the phase transitions at e^ = 0, 
g = Qc and e^ = e^, f7 = found above. By familiar arguments 23| it is expected that 
the universality class of the transition along this entire line will be the same as that of the 
e^ = e^, (? = point, and the e^ = Q, g = Qc point is a singular limit. 
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FIG. 6: Phase diagram of the gauged height model Zu in H2.23() for S = 1/2; the same phase 
diagram also applies to the height model Z'^ in ()2.24|) , or to the exactly equivalent compact QED + 
charge 2 Higgs scalar theory Zcqed in H2.25() . There is no translational symmetry breaking in the 
'deconfined Higgs' phase. The phase with translational symmetry breaking is confining; however, 
it is not completely ruled out that a portion of it could be deconfining, although no such deconfined 
phase with broken translational symmetry appears in Refs |14l. \2nL l2fi|| . Most of the translational 
symmetry breaking is in the columnar pattern of Fig^; however for very large e^ there appears 
to be 26] a change by a first order transition into the plaquette state of Fig^j-this is not shown. 
The entire phase transition line bounding the deconfined, Higgs phase (apart from the single point 
e^ = 0, g = Qc) is described by the fully-frustrated Ising model in a transverse field, Zj in (|2.27|) : 
this transition is expected to be in the D = 2> XY universality class, with the high temperature 
phase of the XY model corresponding to the Higgs phase. The transition at e^ = 0, 5 = Qc is also 
in the D = 2> XY universality class, but it is inverted with respect to the previous one — now the 
low temperature phase of the XY model is the Higgs phase. 

III. DOPED ANTIFERROMAGNETS 



We now discuss recent work[8| which examined the consequences of hole-doping the state 
with broken translational symmetry in FiglHl A theory for hole-doping the deconfined Higgs 
phase of Fig IHl does not exist and would be an interesting direction for future work. 

The study was carried out in a particular large N limit in which the columnar spin- 
Peierls state of Fig^ was a ground state in mean-field theory. Further, the ordering in the 
state was fully developed, and the ground state was well separated from any possible phase 
transition to a deconfined Higgs phase or to a magnetically ordered state. This is probably 
not a realistic situation for the insulator, but it is expected that the results so obtained are 
generic for doping an arbitrary spin-Peierls phase. The results of the study have also been 
reviewed recently [9|, and so we will limit ourselves here to highlighting the main results and 
mentioning some experimental implications. 
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FIG. 7: Ground state phase diagram of a doped antiferromagnet, adapted from Ref [^. A^ is 
the symmetry of spin rotations and it is broken in the hatched region, while C is the translational 
symmetry which broken in the shaded region. At zero doping, the physical system has Neel order 
and so ^A is broken. By using some frustration in the spin Hamiltonian (or considering models 
with a larger spin symmetry group (large N)), we move the system across the quantum critical 
point X into state with columnar spin-Peierls (bond-charge-density) order. The states obtained by 
doping this state are sketched: the thick and dashed lines indicate varying values of the bond charge 
density, while the circles represent site hole density. All the states with non-zero hole density are 
superconducting for large enough A^, but could undergo phase transitions to insulating, anisotropic 
Wigner crystalline states as N is reduced. 



As in Refs p, |7[ , the main effect found was that for reasonable values of the microscopic 
parameters, the holes prefer to segregate in one- dimensional striped structures with a finite 
density per unit length — so as the net hole concentration tends to zero, the stripes move 
farther apart but maintain a finite hole density per unit length per stripe. However, the 
existence of a background of spin-Peierls order makes the details of the hole configurations 
quite different from earlier work p, |7|, \2j\ in experimentally significant ways, as we detail 
below. A schematic for the evolution of the hole configuration is indicated in the phase 
diagram in Fig[7| 

The unit cells of the charge ordered states in Fig [7| have size p x 1; within each unit cell, 
the holes are concentrated in a region of size g x 1. Both p and q were always found to be 
even, with q remaining fixed while p — * cxd as the doping concentration 5 — > 0. Significant 
properties of these even-width stripes are: 

• The holes density per unit length in the g-width region is not unity (as found in earlier 
theories jalll) but various continuously dependent upon microscopic parameters. For 
reasonable choices we can obtain a hole density per unit length of around 0.5 (for g = 2 
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FIG. 8: Charge ordering wavevector (l/p) as a function of the doping, 6 (from Ref [3]) 



this corresponds to a ladder with ^1/4 hole per site), but it is never pinned at exactly 
this value for any finite range of doping, 6 i.e. the stripes are not incompressible. 

There are strong pairing correlations between the holes in each g-width region. This 
leads to an anisotropic superconducting ground state, and should allow good metallic 
conduction above the superconducting transition temperature. These characteristics 
are consistent with observations 28| on La2-x-j/NdySrj:Cu04. 



The evolution of the charge ordering wavevector [K = l/p) as a function of doping, 5, 
is shown in Fig|H| Note that this wavevector is always quantized at values which equal 
l/(even integer). The plateaus are quite small for small 6, but become progressively 
broader as the even integer decreases. In particular, as is observed experimentally, 
a large plateau at wavevector 1/4 is present and this emerges naturally from our 
theory-we regard this as significant. Plateaus at smaller wavevector (say 1/6) have 
not been experimentally observed, but the resolution of current experiments is not 
precise enough to distinguish the staircase- like evolution in Fig |H1 from a continuous 
variation before reaching the plateau at 1/4. The hole density per unit length of stripe 
varies continuously within each plateau but jumps discontinuously as the transition is 
made from one plateau to the next. 

Note that, in FiglHl before reaching the d-wave superconducting state with p = 1, there 
is also a small plateau at p = 2: this is a state with uniform site-charge density, but 
with the bond-charge density of Fig[T^ coexisting with d- wave-like superconductivity. 
Such a state has not been observed, and it would be interesting to look for one. 

All of the stripes found in this computations are bond-centered, at least in the re- 
gion where the symmetry of spin- rotations, A4, remains unbroken; we cannot rule 
out the possibility that there will be a transition to site-centered stripes, like those 
found earlier ja |7|, once Ai is broken. It would be useful for future experiments to 
detect this distinction between bond-centered and site-centered stripes. Let us dis- 
cuss this difference more precisely for the case p = 4. For the site-centered stripe, 
the hole density per unit length in each column of sites takes values pi, p2, Ps, P2 
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before repeating periodically (where pi_3 are some three distinct densities); this is the 
configuration usually assumed in most experimental papers. In contrast for the bond- 
centered state, these densities take the values pi, pi, p2, P2, and this also appears to 
be compatible with existing observations. In the computation here the bond-centering 
is important for enhancing pairing correlations, which are responsible for supercon- 
ducting/metallic transport in the direction parallel to the stripes. We note that a very 
recent NQR experiment i29] indicates a charge distribution which has some features 
consistent with the bond-centered state: they find only two inequivalent sites associ- 
ated with the stripes, and a density per site in the hole-rich region which is only about 
0.18-0.19 (the fully segregated site-centered stripe would have pi ~ 0.5, P2 = Pa = 0, 
while the fully segregated bond-centered stripe would have pi ~ 0.25, p2 = 0; clearly 
the bond-centered value is more compatible with observations). 
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